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In soliton models expressed in terms of the non-linear chiral field, the electric current has an 
anomalous gauge-field contribution as the baryon current does. We study the spin polarized 
Skyrmions coupled with the electromagnetic field via the gauged Wess-Zumino term and calcu- 
late configurations of the Skyrmion and the gauge field with boundary conditions to ensure the 
physical charge number for baryons. 

Although the electromagnetic field via the gauged Wess-Zumino term affects physical quantities 
f~^ , in small amounts, we find that the magnetic field forms a dipole structure owing to a circular 

f^ . electric current around the spin quantization axis of the soliton. This is understood on an analogy 

^\l ' with the Meissner effect in the super conductor. The electric charge distributions turn out to have 

characteristic structures depending on the total charge, which suggests the intrinsic deformation of 
baryons due to orbital motions of the constituents. 
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^ , INTRODUCTION 

> 

^' . n . . 

f~^ Of recent years the proton spin puzzle has attracted much attention. As reported in Ref. ^1], the quark spin is 

^~^ found to be responsible for only a small amount of the proton spin. The gluon spin contribution is also estimated 

^^ from the global analysis of the parton distribution function _2J, but has yet large uncertainties. The experiment of 
^~l ■ polarized proton collision at RHIC in progress is expected to give some clue to the gluon spin content. Such being the 
f— ^ r situation, it is desirable to study the possibility that the orbital angular momentum may contribute to some extent to 
the proton spin and the orbital motion may induce intrinsic deformations of the baryons. Non-spherical components 
in baryons have been discussed in the context of the color magnetic interaction^, the electromagnetic transition^ 
and the generalized parton distribution 5] . 
^-1, Deformed baryon states are investigated also in soliton models[6|, [tI [a]. Since the topological soliton is quantized 

to the fixed (iso-)spin states by rotation, the intrinsic deformations directly correlate with orbital motions in this 
description. Considering that the shape of the charge density is probed via the photon, we study the Skyrme model 
coupled with the electromagnetic field. Note that there are two kinds of terms which are brought about by the 
coupling of the gauge field; non-anomalous terms through covariant derivatives and anomalous terms through the 
^ ' gauged Wess-Zumino (WZ) term. 

The gauge fields are minimally incorporated through the covariant derivatives. Using variation of an action obtained 
by the minimal replacement, the authors of Ref. |9( computed configurations of the gauge field as well as the soliton 
and estimated the magnetic moment. They found that the dipole magnetic field is generated around the soliton 
and the magnetic moment has reasonable value. However the topological current is not taken into account in their 
variational equations. 

The topological baryon current is naturally incorporated in the Maxwell equation if one considers the gauged WZ 
term. The gauge fields affect the system not only through the covariant derivatives but also through the gauged WZ 
term, which is not obtained by the minimal replacement of derivatives in the WZ term. The gauged WZ term is 
designed to account for the non-conservation of the axial current by the anomaly and to describe correctly anomalous 
processes like 7r° -^ 27 in the pion sector. This term possibly influences the system also in the soliton sector. Above 
all, the gauged WZ term is indispensable to assure the gauge invariance of conserved currents like the baryon current 
and the electric current as well. Actually, the gauged WZ term provides an anomalous contribution lCr| of the magnetic 
field to the isospin charge of the soliton. 

Although the electromagnetic field is usually treated perturbatively owing to the small coupling constant, it is 
worth estimating the anomalous contribution of the gauge field as the WZ term is essential for the anomalous baryon 
current. It is possible that the gauged WZ term exerts significant infiuences on the physical quantities of the isospin 
fixed solitons. 
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In this work, we adopt an action including the gauged WZ term and construct spin polarized Skyrmions with 
a proper electric charge imposing boundary conditions upon the chiral field and the gauge field. We calculate the 
configuration of these fields and study the properties of the spin polarized Skyrmions, particularly characteristic 
spatial structure of the electric charge, current and the magnetic field. Through these results, we discuss the intrinsic 
deformations of the charge density as well as the significance of the effects of the electromagnetic fields from the 
gauged WZ term in the soliton model. 

THE SKYRME MODEL WITH THE GAUGED WESS-ZUMINO TERM 

We present a model written in terms of the chiral field and the gauge fields. To make it manifest how the gauge 
fields are incorporated in the system, we consider a chirally gauged model first in a general context and then restrict 
ourselves to a model coupled with the electromagnetic field only. The action of the chiral field consists of the non- linear 
sigma term, the Skyrme term with pion mass term and the Wess-Zumino term as 
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The Wess-Zumino term, /wz, is expressed compactly by the use of differential forms as 

*^' ^ Tr ((uUufAi^ - (uUU){Ai^dAi^ + dAi^Ai^ + Al) - l-{uUUAi^f + {uUu)^U^AnUAi 



487r2 
- dAi^{uUu)U^AnU - {uUu)A^U^ArUA^ -{U^AnU){Ai^dA^ + dAi^A^ + Al) - ^{U^AnUAi^f - (p.c.) 

with the abbreviation p.c. for parity conjugate: U ^^ U'^,Ai^ «-> Ar. The covariant derivative for the chiral field is 
defined as dx'^Df.U = dU + ArU - UAi^. 

This action except the Wess-Zumino term and the pion mass term has the gauge symmetry and is invariant under 
the gauge transformation: 

U -^ 9RUg^\ Al ^ 5l(Al + d)g^\ Ar -^ gniA^ + d)g^\ (1) 

The Wess-Zumino term is constructed so as to reproduce the chiral anomaly attended by this transformation and 
hence cannot be obtained by the minimal replacement of differential operators which makes an action gauge-invariant. 
The Noether currentjlflj associated with the vector symmetry, i.e. ^l = 3r, generated by Q is 

jM = _ ^Tr iQiUD'^U'' + (p.c)) + — ^Tr ([iQ, UD^U^][UD^'U\UD''U''] + (p.c.)) 
8 8e"^ ^ ' 

- i^^^^^'^Tr ^Q[iUDU^)l,^ - {UD^U^F^x, + ^UF^^pU^ - (p.c.)], (2) 

where F^, Fi, are the field strengths for the gauge fields Ar, Al respectively. If we set Q = 1/A^c to get the baryon 
current, only the last term — the contribution from the WZ term — survives and we find the baryon number is 
equivalent to the winding number. Therefore we can identify a baryon with a topological soliton in this model. 

Now we restrict ourselves to study of the SU(2) chiral field and the Uem(1) gauge field, Al = Ar = ieQEuAfj.dx'^ . 
The electric charge matrix satisfies the Gell-Mann-Nishijima relation: Qem = {t^ + l/-^c)/2. A soliton of a unit 
baryon number is realized in the hedgehog form : U = exp(zT • rF{r)) with the boundary conditions F{r = 0) == 
TT, F{r = oo) = 0. The electric potential is assumed to be spherical 

eA° = Vir) 

for simplicity, and the angular dependence of the vector potential is fixed as 

eA = h(r) sin^ BVcj) 

to match the Ampere law. Here, 9 and (j) are the polar and azimuthal angle respectively. All these do not depend 
on time and such a static configuration is always taken by the gauge transformation Q. Even a (time-dependent) 



collective rotation of the chiral field can be described by taking g^ = gi, as the rotation matrix, if one lets eA° take 
a constant value in proportion to the angular velocity. 

The total Lagrangian density including the kinetic term of the electromagnetic field is written in terms of F{r), V{r) 
and h{r) as 
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where Sp = sin F{r) and Cp — cos F{r) . After the angular integration and scaling r, V, 7X1^^ by eF^^ to be dimensionless, 
we are led to the variational equations; 
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the Gauss law V ■ E — j'^-^ with the angular averaged charge density, 
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and an equation deduced from the Ampere law V x B = Jemi 
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with the fine structure constant a = e^ /An. 

Note that the magnetic field B (or h in r.h.s of Eq. Q) is instrumental in producing partly the electric charge 
Jem ^^ ^ source of the electric field E and that E contributes the electric current j'^j^ as a source of B. This 
complementarity comes from the gauged WZ term in the action which provides the anomalous couplings between 
electromagnetic field and the chiral field like E ■ B tt'^/F.^. The term brings to j^^^ an anomalous current of the dual 
field strength as e^'^^POxApd^ir^ /F^ with the non-vanishing pion field d^Tr°/F^ ~ i Tr T^{Ud^W - Wd^U) in the 
soliton sector. In addition to this, here we emphasize that the topological term of the baryon density is correctly 
incorporated in the charge density due to the gauged WZ term, as one can see the first term in r. h. s. of Eq. Q. 

We solve the coupled equations lEl)-© imposing boundary conditions for the fields. 



F(0) = TT, 

l/'(0) = 0, 
h{r -^ 0) oc r^, 
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where Z is the charge number of the soliton (to be set at 1 or 0) and fi^ — rn^ — 2V^/3 is a modified pion 
mass|3, lllf caused by 'rotation' of the soliton. The angular velocity is turned to an asymptotic constant Voo by a 
gauge transformation as mentioned previously and thus non-zero Voo corresponds to a 'rotation' of the soliton. In 
practical computing, we impose a boundary condition on V'{r) instead of V{r) without any tuning of Voo and confirm 
finite values, actually. 

However, the rotation matrix does not cover the whole SU(2) group because the generator Qem of Uem(1) gauge 
group includes only the z component of the isospin matrix. Accordingly, A cannot be excluded from nucleon states 
without the Casimir projection as performed in Ref. ,13] . In this sense, our treatment is not equivalent to the 
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FIG. 1: Electromagnetic fields with the parameter set II are plotted as a function of radial distance r. Bold lines are for 
the charged soliton and thin lines are for the neutral one. The electric potentials (normalized by efV) with the asymptotic 
constants subtracted are shown by solid lines and — /i(r)'s which determine the magnetic field are by dashed lines. 

collective quantization of rotation fully, but to one-dimensional cranking. Nonetheless, the z component of isospin 
(as the electric charge minus one half of the baryon number) is assured of a quantized value ±1/2 for the charged 
and neutral soliton (Z = 1, 0) only by the boundary conditions for F and V , though this does not mean a topological 
quantization of isospin in contrast to that of the baryon number. Since the asymptotic value Voo affects the boundary 
condition of i^(r), the self-consistent approach we took to solve the differential equations corresponds to the variation 
after projection. Furthermore, the z-component of the spin takes Tl/2 for the charged and neutral solitons, since the 
spatial rotation means the inverse iso-rotation for the hedgehog configuration and the photon contribution to the spin 
is negligible which we confirm numerically. Therefore we regard the soliton as a spin polarized nucleon with possible 
admixture of the A component. 



FIELD CONFIGURATIONS AND CHARGE DISTRIBUTION 



We solve the differential equations with the boundary conditions for proton [Z — 1) and neutron {Z = 0) using the 
relaxation method. Three sets of parameters -FV and e are chosen while m.^- is fixed at the physical pion mass. Set 
I: (e = 5.37, Ft, = 185 MeV) is determined so that the charge radius of the proton is reproduced with physical pion 
decay constant. Set II: (e = 5.45, F.„ = 129 MeV) and Set III: (e = 4.84, F^ = 108 MeV) are taken from Refs. ^ and 
[l3], respectively. We find that for Set I, tq such that F{ro) = tt/2 decreases by a few % compared with that of Set II. 
While this means F{r) for Set I is squeezed, F{r) for Set III, on the contrary, swells compared to Set II. Although the 
boundary condition constrains the fields to have the different asymptotic behavior, profiles of F{r) with Sets II and 
III are almost the same as that obtained in Refs. 12i| and 13] respectively. The configurations of the electromagnetic 
field are also determined simultaneously, as shown in FIG. ^ The electric potentials display correctly the Coulombic 
behavior for large r, and the smeared peak near the origin means finite structures of soliton's charge unlike a point 
particle. The function h{r) gives the magnetic field as 
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and the radial distance of their peak determines the size of the dipole magnetic field. The spatial distribution of the 
fields discussed in detail later. Physical quantities with these parameter sets are shown in TABLE Q] 

For Sets I and III, the mass of solitons deviates from the experimental values, but these receive corrections of meson 
loop and higher chiral-order terms |l4|. Mass differences between proton and neutron are of proper magnitude for all 
sets of parameters but show the opposite sign. This is because the spherical ansats for the electric field is oversimple 
to give the correct sign, though we take the electric field into account even for the neutron. In fact, the energy of the 
electric field contributes about 0.2% to the proton mass which is reasonable for the mass difference, while it almost 
vanishes for the neutron case. 



This setting is also reflected in the charge radius defined as {r'^)ch — J d^x r^ Jem/^ with the charge density, 



Jem 



4:'K^r^ 
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_^ri + 4F- + 4^)5^ysin^,, (8) 



which is obtained with assignment Q = Qem in Eq. ^ and scaling to the dimensionless variables. For the proton, the 
positive charge density is piled up and the charge radius is well reproduced comparatively. For the neutron, however, 
the positive and negative charge densities compete with each other and the charge radius of neutron is overestimated 
several times as much. 

The main contribution to the charge density, ((HJ, comes from the first term of the baryon- number density near the 
origin and from the last term of the isospin charge for the medium range. Since the latter has a remarkable polar-angle 
dependence, the spatial distribution of the electric-charge density presents a deformed or a toric structure with a core 
according to the isospin, as shown in FIG.|21 

The baryon number density depends also on the polar angle through the gauge field. 
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but it is almost spherical because the small fine structure constant suppresses the magnitude of field h. On the other 
hand, the isospin charge is mainly contributed from the non-linear sigma term proportional to SpV sin 9 as is seen 
in Eq. ©. This term affects the electric charge with the opposite sign for charged and neutral solitons, especially 
in the region such that F{r) ^ 7r/2 and svaO ~ 1. As seen in FIG.|31 the solid curves of the charged soliton change 
according to the polar angle and smoothly merge into the dashed lines of the neutral one at the polar region. Although 
the angular-averaged density of the electric charge agrees with that obtained in Ref. |lj|, the 9 dependence makes 
an oblate shape for the charged soliton and a prolate positive-charge surrounded by a toric negative-charge for the 
neutral soliton. The deformed charge density is interpreted as a consequence of the centrifugal force accompanying 
the rotation of soliton. The radius of the negative-charge torus is about 0.7 fm with the parameter Set II, which 
is governed by the distance tq such that F{rQ) ~ 7r/2. Since such a deformed structure cannot be realized for the 
nucleon with spin 1/2, the density distribution is interpreted to stand for an intrinsic deformation of the nuclcon and 
also to refiect the admixed A component. 



The magnetic moment /i is evaluated from the electric current as 2 / d'^x r x Jem with 
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which is equivalent to the estimation from asymptotic behavior of the magnetic field, linir^oo rh{r)/a, owing to the 
Ampere law. Our results in unit of the nuclear magneton are comparable to the experimental values including their 
sign. 

The reasonable estimation of the magnetic moment results from the dipole magnetic field generated by the current 
of soliton. Actually, the magnetic field expressed by Eq. has a dipole structure induced by the circular current Eq. 
©, as shown in FIG. 01 These field configurations are the same for both of the charged and neutral cases, while the 







mass 

[MeV] 


[fm^] 


[fm] 


[fm] 


[Miv] 


PA 


QirNN 


a 

[MeV] 


Set I 


P 


1326 


0.757 


0.416 


1.128 


3.27 


0.616 


6.08 


31.4 




n 


1325 


-0.556 


0.415 


1.164 


-2.84 


0.616 


6.25 


31.3 


Set II 


V 


933 


0.773 


0.563 


1.118 


3.95 


0.572 


8.20 


34.1 




n 


931 


-0.446 


0.563 


1.136 


-3.28 


0.571 


8.35 


34.0 


Set III 


P 


889 


0.890 


0.708 


1.086 


5.53 


0.672 


11.9 


42.8 




n 


887 


-0.385 


0.707 


1.124 


-4.79 


0.671 


12.1 


42.7 


exp 


P 

n 


938 
940 


0.757 
-0.116 


0.801 


0.81 


2.79 
-1.91 


1.26 


13.5 


~30 



TABLE I: Static properties of the nucleon described by the soliton with several sets of parameters. {r^)ch, {i''^)b and {r^)M are 
the mean square charge, baryon number and magnetic radii. 



spin of the neutral (charged) soliton is polarized (anti-)parallel to the z axis. The radius of the circular current is also 
controlled by the scale ro satisfying F(ro) ^ 7r/2. The spatial distribution of the magnetic field is understood also 
on an analogy with the Meissner effect in the super conductor. The coefficient of the vector potential in the current, 

( ^ — ^{(<9pi^)^ + sin^ F(9^0)^} 1 sin^ Fsin^ 0, plays a role of a square 'mass' of the gauge field, because this factor 

can be read as the coefficient of (e^^)^ in the action. In case of the hedgehog configuration, the large 'mass' region 
forms a torus whose core, as one can see from the factor sin i^sin 0, is located at the circle around the z-axis such 
that F(r) = 7r/2 and d = 7r/2. The configuration of the magnetic field is determined so as to avoid the large 'Meissner 
mass' region to save the total energy. As a consequence of this fact, the magnetic field coils around a genus of the 
torus. 

The axial coupling tends to be one half of experiments, which agrees with Refs. |l2L ll.lj. The pion-nucleon sigma 
term evaluated from the symmetry breaking term results in the adequate value as shown in TABLE Q] 

In any case, contributions of the electromagnetic field via the gauged Wess-Zumino term are controlled by the 
fine structure constant and hence amount to 1 % of the physical quantities, at most. Therefore the complementary 
dynamics of the electromagnetic field is not conspicuous in the present situation. However, this anomalous contribution 
might be influential in high dense matter. 

CONCLUDING REMARKS 

We have studied the effects of the electromagnetic field coupling to the soliton configuration via the gauged Wess- 
Zumino term. We pointed out that the gauged Wess-Zumino term provides the Maxwell equation with anomalous 
terms in the presence of the soliton. Unless this term is considered, the topological baryon current cannot correctly be 
incorporated in the charge density of variational equations. Furthermore non- vanishing pion fields make the anomalous 
coupling -K^E ■ B effective in the soliton sector, and the electric and magnetic fields contribute complementarily to 
the source current of each other. Such a current from the dual field strength is interesting theoretically, but in reality 
these effects are suppressed by the small coupling constant. 

Because the anomalous contribution of the electric field is small, the spatial distribution of the electric current is 
mainly determined by the non-anomalous term of the isospin current. From this, it is concluded that the electric 
current has the toric configuration, and accordingly that the magnetic field has a poloidal structure coiling around 
the circular current according to the Ampere law. This dipole magnetic field is understood as a consequence of the 
Meissner-like effect. Actually, the Meissner 'mass' of the gauge field is estimated from the chiral angle F as sin^ F sin^ 6 
corresponding to the condensate of the Cooper pairs. Contrast to the superconductor, the large mass region does not 
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FIG. 2: The charge distribution of the charged soliton in the xz plane (left) and the neutral soliton in the xyz space (right) 
calculated with the parameter Set II. Note that the soliton spin is polarized along the z axis. Gray and black dots are positive 
and negative charges. The charge density of the charged soliton is deformed to an oblate shape and that of the neutral soliton 
has a Saturnian structure. 
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FIG. 3: The charge densities as functions of the radial 
distance r with several polar angles 9. 
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FIG. 4: The arrows in the xz plane show the magnetic 
field, which is axial symmetric. Dots in the negative x 
region represent the current density running around the 
z axis clockwise. 



extend uniformly but forms a finite torus, and hence the magnetic flux quantization cannot be achieved in our case. 

As discussed above, the symmetric axis of the torus of the field configuration is to be interpreted physically as a 
spin-quantization axis of the soliton. This leads us to a reasonable consequence that the magnetic moments evaluated 
from the asymptotic form of the poloidal magnetic fields are (anti-) aligned with the spin of the charged (neutral) 
soliton. The theoretical value of the magnetic moment derived from the variational solution is of the same order as 
the experimental value of the nucleon. 

We have found that the electric charge distributions are distinctive for their polar angle dependence. The equatorial 
swelling of the isospin charge makes the difference while the spherical core of baryon number density is common to 
charged and neutral solitons. Our results show that the charged soliton has the oblate shape of the charge density and 
the neutral soliton has the toric negative charge surrounding the prolate shape of positive charge. These distributions 
are consistent with the pion cloud rounding the nucleon or may imply a u-d diquark core accompanied by the other 
valence quark. In any case, the characteristic distributions of the charge density suggest the intrinsic deformation of 
baryons and the isospin-dependent quadrupole moment, which might be observed in the A^A transition. 

As in the dynamo theory for the planetary magnetic field, a poloidal magnetic field may induce a toroidal magnetic 
field through the rotation of the soliton, and vice versa, even though the physical scale is quite different. It is our 
future work to take account of this effect by extending the variational space of the magnetic field. The non-radial 
electric field is also important to estimate the mass of the proton and the neutron and to discuss a correlation between 
the orbital motion and the intrinsic deformation. 
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